The role of orbital dynamics in spin relaxation 
and weak antilocalization in quantum dots 
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We develop a semiclassical theory for spin-dependent quantum transport to describe weak 
(anti)localization in quantum dots with spin-orbit coupling. This allows us to distinguish different 
types of spin relaxation in systems with chaotic, regular, and diffusive orbital classical dynamics. We 
find, in particular, that for typical Rashba spin-orbit coupling strengths, integrable ballistic systems 
can exhibit weak localization, while corresponding chaotic systems show weak antilocalization. We 
further calculate the magnetoconductance and analyze how the weak antilocalization is suppressed 
with decreasing quantum dot size and increasing additional in-plane magnetic field. 
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Weak localization (WL) and antilocalization (AL) are 
classic examples for quantum interference and spin-orbit 
(SO) interaction effects on the conductance in low- 
dimensional electronic systems Q, Q . Very recently, par- 
ticularly weak AL has been reconsidered in a number of 
experiments since AL can be employed as a probe of SO- 
induced spin dynamics and relaxation phenomena. Mea- 
surements have been performed both for GaAs- and InAs- 
based two-dimensional(2D) electron gases Q, as well as 
for ballistic bismuth and GaAs [5j cavities. While 
SO scattering in extended disordered systems is well un- 
derstood 0, the latter experiments address the timely 
question of how quantum confinement of the orbital mo- 
tion affects spin relaxation in clean ballistic quantum dots 
where the elastic mean-free path is much larger than the 
system size. Considerable related progress has also been 
made theoretically in treating spin relaxation and the in- 
terplay between SO and Zeeman coupling in quantum 
dots [a] including random-matrix theory (RMT) 0, Q • 
However, the RMT results apply only to chaotic quan- 
tum dots and contain geometric parameters which must 
be obtained by other means for a given system. 

Here we present an alternative, semiclassical theory 
for the spin-dependent magnetoconductance of quantum 
dots, i.e. a semiclassical Landauer formula including 
spin, and apply it to describe weak AL in 2D confined 
systems. This approach allows to uncover the inter- 
relation between orbital dynamics and spin evolution in 
a transparent way, and it is rather generally applicable to 
quantum dots with different type of classical dynamics, 
e.g., chaotic and regular. Remarkably, we find signifi- 
cant qualitative differences in the spin relaxation times of 
chaotic, integrable, and open diffusive systems: Spin re- 
laxation for confined chaotic systems is much slower than 
for diffusive motion; moreover, for a number of integrable 
geometries we even find a saturation, i.e., a certain spin 
polarization is preserved. Furthermore, we examine the 
effect of the system size and of an additional in-plane 
magnetic field on the resulting AL. 

Our study is based on the semiclassical Landauer for- 



mula 0, that we generalize to systems with SO and 
Zeeman interaction. To this end, we extend techniques 
for spin semiclassics, recently developed for the density 
of states 0,0], to quantum transport. We consider a 
Hamiltonian linear in the spin operator s, 



H = H (q, p) + /is-C(q, p), 



(1) 



where C(q, p) is a vector function of the position and 
momentum operators q, p, which may include an ex- 
ternal (inhomogeneous) magnetic field. For a large 
number of systems of interest, and usually in experi- 
ments, fts|C(q,p)| <C H , even if the spin-precession 
length is of the order of the system size. Here s is 
the particle spin, and the phase-space functions with- 
out a hat denote the classical counterparts (Wigner- 
Weyl symbols) of the respective operators. As a con- 
sequence of the above inequality the influence of spin on 
the orbital motion can be neglected. Thus Ho deter- 
mines the classical trajectories 7 = (q(i),p(i)) which, 
in turn, generate an effective time-dependent magnetic 
field C 7 (t) = C(q(i), p(i)) acting on spin via the Hamil- 
tonian Hy (t) = H s • C 7 (t) . Hence the spin dynamics can 
be treated quantum-mechanically in terms of a (time- 
ordered) propagator K 7 (t) =T exp[— i J Q dt's ■ C 7 (£')]. 

To derive a semiclassical expression for the spin-depen- 
dent conductance of a quantum dot, we start from the 
Landauer formula in two dimensions relating the two- 
terminal conductance G = (e 2 /h)T to its transmission 
coefficient 0] 



N' N s 

r = EE E 1*- 

n=l m—1 a.a' — — s 



(2) 



The leads support N and N' open orbital channels m 
and n, respectively, and we distinguish 2s + 1 spin po- 
larizations in the leads, labeled by a — —s, . . . , s. In J2J, 
tntr'.mtr is the transition amplitude between the incom- 
ing channel \m,a) and outgoing channel \n, a'); a corre- 
sponding equation holds for the reflection coefficient 1Z 
satisfying the normalization condition T+1Z = (2s+l)N. 
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A semiclassical evaluation of the transition amplitudes, 
starting from a path-integral representation of the Green 
function, yields 01 ( see m f° r the spinless case) 



tna' ,m<7 — ^ ^ (-^"7)0"' cr A-y exp [ ^^7 

7(n,m) 



(3) 



given as coherent summation over classical paths at fixed 
energy corresponding results hold for the reflection 
amplitudes in terms of back-reflected paths. The sum 
runs over classical trajectories j(fi — ±n, fh — ±m) that 
enter (exit) the cavity at "quantized" angles 6^ (On)- 
For hard- wall boundary conditions, sinOm = mn/kw 
and sinO fi = fnt/kw', where k is the wavenumber, and 
w, w' are the lead widths. In J3J), S 1 (=hkL 1 for billiards) 
is the action along 7 with time T 7 and classical weight 



The entire spin effect is contained in the matrix 



elements (-ftT 7 ) (T ' cr of the spin propagator K 7 = if 7 (T 7 ) 
between the initial and final spin states. 

Inserting © into (J2J) we derive the semiclassical Lan- 
dauer formula for spin-dependent magnetotransport (in- 
cluding SO and Zeeman interaction): 



T 



E E 

nm 7(n,rn) 7'(n,m) 



* .(i/R)(5 7 -Sy) 



(4) 



The orbital contribution of each pair of paths is weighted 
by the spin modulation factor 



M-y, Y =Tr(X 7 £t,). 



(5) 



where the trace is taken in spin space. 

WL and AL effects are obtained after energy average 
of T(E, B) for ballistic quantum dots (subject to an ex- 
ternal arbitrarily directed magnetic field B). The leading 
contributions after averaging Q for a chaotic cavity with 
time-reversal symmetry, i.e., B = 0, are as follows: 

(i) The classical part consists of the terms 7' = 7 , 
for which the rapidly varying energy-dependent phase in 
the exponent of (0} disappears. Then the modulation 
factor is M ln = Tr (K 7 K^) = 2s + 1, independent of SO 
interaction, and reduces to the trivial spin degeneracy. 

(ii) The diagonal quantum correction is defined for the 
reflection only. It contains the terms with n = m and 
7' = 7 _1 , where 7 -1 is the time-reversal of 7 0. Again, 
the orbital phases of the trajectory pair cancel; however, 
the modulation factor is .M 7i7 -i =Tr [K~ ). 

(iii) The loop contribution comes from pairs of long or- 
bits that stay close to each other in configuration space, 
thereby have nearly equal actions, and hence persist upon 
energy average. One orbit of the pair has a self-crossing 
with a small crossing angle, thus forming a loop, while 
its partner exhibits an "anticrossing" . Outside the cross- 
ing region the orbits are located exponentially close to 
each other: the paths are related by time reversal along 
the loop and coincide along the rest of the trajecto- 
ries [lfj, [l6( . We have computed the modulation factor 



for 7 and 7' neglecting the crossing region and found 
M 7i7 < — Tr {K 2 ), where / is the loop segment of 7 0|. 

For spinless particles in chaotic quantum dots the three 
contributions to the averaged transmission and reflection 
yield 0, for N = N' > 1, (i) Tj 0) = = N/2, 
(ii) SK { d l g = 1/2, and (iii) ST^ p = 5Tl^ op = -1/4, in 
agreement with RMT. Here the superscript refers to zero 
spin and zero magnetic field. 

In the following, we consider the case of an additional 
uniform, arbitrarily directed magnetic field B in the pres- 
ence of SO interaction. Besides the Zeeman interaction, 
the field component B z perpendicular to the cavity gen- 
erates an additional Aharonov-Bohm (AB) phase factor 
ip = cxp(i47rA 7 _B z /$o) i n the diagonal and loop terms in 
Eq. Here, A 7 = J A • dl/B z is the effective enclosed 
area for the diagonal (loop) contributions accumulated 
along the orbits 7 (loops I) neglecting bending of the or- 
bits, and <!>o = hc/e is the flux quantum. 

For broken time-reversal symmetry, e.g., by the per- 
pendicular B-field, M. ln i should be calculated directly 
from ©. We then introduce a generalized modulation 
factor, M v = M. lri iLp, which is distributed according to 
a function P(A4 ip ; L, B), where L is the trajectory (loop) 
length in the diagonal (loop) contribution, and the B- 
dependence includes both the AB phase and the Zeeman 
interaction. Thus we can define a spin modulation factor 
Ai v (L; B) averaged over an ensemble of trajectories with 
fixed length L. 

In chaotic systems, the length distribution is given by 
exp(— L/L csc ) 0], if the escape length i osc = ttPl c /(w + 
w'), the average length the particle traverses before leav- 
ing the cavity of area A c , is much larger than Lb, the 
average distance between two consecutive bounces at the 
boundaries. It can be shown [ToL ll4T | that the relevant 
distribution of loop lengths is determine d by the same 
exponent. As in the case without spin [9|,ll0(, the prod- 
uct of the AB phase and spin modulation factors in Eq. 
(0} can be eventually substituted by its average (Ai v ) l 
over L and pulled out of the sum. Thereby we obtain, 
as relative quantum corrections for the spin- and B-field- 
dependent transmission and reflection, 



^ diag /^^ g = «w<«4 op 



6T loop /6T^ p 



1 fOO 

= (M V (B)) L = — dL e - L ' L - M V (L; B) . (6) 

^csc Jo 

Note that current conservation, i.e., STZ^g + <57^ioop = 
— 6T\ oop , is fulfilled in the semiclassical limit N,N' ^> 1. 
In the absence of SO interaction, we have A^ v (i;B) = 
(2s + 1) cxp(-S 2 L/ J L b ), where B = 2V27rB 2 A /$ and 
A is the typical effective area enclosed, and the usual 
Lorentzian S-profile [j|E3 is recovered by Eq. JHJ. 

In the case of SO interaction, the quantum correc- 
tions depend on the modulation factor .M V (L;B), 
which characterizes the average spin evolution of a tra- 
jectory ensemble and can be easily determined from clas- 
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sical numerical simulations. Without B-field, A4(L) = 
M. ip (L\ 0) changes from A^( 0)~2s + 1 to the asymptotic 
value M.(oo) = (— l) 2s I4j; i.e., for s = 1/2 an initial 
polarization (Ai(Q) — 2) becomes completely randomized 
(M(oo) = —1) owing to SO interaction, if the particle 
motion is irregular (see below). Thus, if the particle 
quickly leaves the cavity (large openings, small L csc ) or 
the SO interaction is too weak, there is not enough time 
for the modulation factor to deviate from 2s + 1, giving 
rise to standard WL. In the opposite limit (large L osc 
or relatively strong SO coupling) A4(L) quickly reaches 
its asymptotic value and, in view of ©, (-M v (0))l — 
(— l) 2s . Hence for half-integer spin the conductance cor- 
rection becomes positive due to SO interaction. This 
phenomenon of weak antilocalization does not exist for 
integer spin. For B ^ we find A^oojB) = 0: 
Both a magnetic flux (destroying constructive interfer- 
ence of the orbital phases) and the Zeeman interaction 
(affecting the spin phases) inhibit AL. 

For a quantitative treatment we must specify the form 
of the SO interaction. In the following numerical analysis 
we consider the spin s — 1/2-case for different quantum 
dot geometries and Rashba SO coupling 0, relevant 
for 2D semiconductor heterostructures. It is described 
by an effective magnetic field C = (2 ctR m e /H 2 ) v x z, 
where ctR is the Rashba constant, m e is the effective 
mass, and v the (Fermi) velocity. In a billiard with fixed 
kinetic energy, C is constant by magnitude and its di- 
rection changes only at the boundary. It is convenient 
to characterize the SO interaction strength by the mean 
spin- precession angle per bounce, Or = 2irL b / Lr, where 
Lr = 2it\v\/C is the Rashba length. 

In the inset of Fig. ^ we plot the modulation fac- 
tor A4(L) for three SO strengths for a chaotic, desym- 
metrized Sinai (DS) billiard (Fig. geometry 2), i.e., 
a prototype of a geometry with hyperbolic classical dy- 
namics. The average was performed over an ensemble of 
10 5 (non-closed) trajectories (in the closed system) with 
random initial velocity directions and positions at the 
boundary. As 9r increases, M.{L) reaches its asymptotic 
value —1 faster. In the main panel of Fig.^we compare, 
for fixed 9r/2it = 0.2, A4(L) for four systems represent- 
ing three different types of orbital motion: 

Chaotic systems: Two representative geometries, the 
DS billiard (curve 2) and the desymmetrized diamond 
(DD) billiard [l8j (curve 3), show up to deviations at 
small lengths nearly the same decay behavior of A4(L), 
indicating universality features for chaotic dynamics. 

Integrable systems: Although Eq. © is valid only for 
chaotic cavities, the average modulation factor M.(L) is 
well defined for other types of motion. Remarkably, we 
find that for the integrable quarter-circle (QC) billiard 
(curve 1) M(L) oscillates (with frequency independent 
of 9r) around a constant saturation value well above —1. 
A systematic analysis shows |l4j that the saturation value 
in the integrable case is system-dependent and decreases 
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FIG. 1: Average modulation factor M(L) as a function of or- 
bit length L in units of bounce length Lb for the quarter-circle 
billiard (curve 1), the desymmetrized Sinai billiard (curve 2), 
the desymmetrized diamond billiard (curve 3), and for an un- 
bounded diffusive system with mean-free path equal to Lb 
(analytical curve 4). The relative strength of spin-orbit inter- 
action is #r/27t = 0.2. Inset: A4(L) for the desymmetrized 
Sinai billiard at different values of On. 

down to —I, indicating spin relaxation, with increasing 
Or. 

Diffusive systems: Unbounded diffusive motion (curve 
4) exhibits fast exponential relaxation; i.e., A4(L) ~ 
3exp[-(6» 2 /3)(L/L 6 )]-l 0, where L b is identified with 
the scattering mean-free path (cf. Eq. (10.12) of Ref. 0]). 

Note that the curves 1-4 almost coincide for L < L b , 
because up to the first scattering event the particle moves 
along a straight line, and different types of dynamics can- 
not be distinguished. On larger length scales, we find 
significant qualitative and quantitative differences in the 
spin evolution in chaotic, integrable, and diffusive sys- 
tems. In particular, the relaxation is strongly suppressed 
for a confined, even chaotic, motion as compared to an 
unbounded diffusive motion with the same Or. This re- 
sult is supported by the following argument: In the limit 
6r ^ 1 the spin movements on the Bloch sphere "mimic" 
the orbital motion to order 0^; i.e., they are bounded 
for a spatially confined system. If higher-order correc- 
tions were neglected, the spin relaxation would saturate 
at L ~ Lb- The further decrease of A4(L), of order 
(L/Lb) Or, is due to a Berry phase acquired by the spin 
wave function 0, • Its effect is similar to that of the 
AB phase. Hence, in a chaotic system without Zeeman 
interaction one finds 

JaT^L; B) ~ e ~CB+~6lfL/L b + e -(£-?|)=L/L b) (?) 

with 9ft = (k /L b 2 ) R l\f2. The further relaxation is due 
to terms of order (L/Lf,) Or ^ eventually renders 

Ai(L) negative and causes AL. For stronger interaction, 
6r ~ 1, the three mechanisms (initial relaxation, Berry 
phase, and further relaxation) work simultaneously and 
cannot be separated (e.g., curves 2 and 3 in Fig. QJ. 
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FIG. 2: Relative quantum correction to the reflection 
SIZ/SIZ^ vs. spin-orbit interaction Or for B = in the de- 
symmetrized Sinai (solid), diamond (dashed), and quarter- 
circle (dashed with circles) billiards with P c /(w + w) = 90. 
Lower left inset: 81Z/81Z^ vs. Zeeman interaction 9z for the 
desymmetrized Sinai billiard. The in-plane field is directed 
parallel (solid) and perpendicular (dashed) to the long side. 
Upper right inset: 51Z/S1Z^ vs. perpendicular magnetic flux 
for the same billiard with 9z = 0. 



where one can switch from WL to AL by tuning the clas- 
sical dynamics from integrable to chaotic. 

The Zeeman interaction, measured by a precession an- 
gle 9z per bounce (analogous to Or), suppresses AL 
(lower left inset). Note the anisotropy in the (in- 
plane) field direction. The upper right inset depicts the 
magnetic-flux dependence. The characteristic double- 
peak structure follows from Eqs. (jfJJl and Q). 

The present semiclassical approach has a wider range 
of applicability, including ballistic integrable systems and 
SO strengths up to Or ~ 1, compared to RMT 0, H, 
which assumes 0r<c1 in the ballistic regime (Eq. (23) of 
Moreover, the RMT results contain free geometric 
parameters that have to be computed separately. 

A corresponding analysis of ballistic conductance fluc- 
tuations [3] with spin appears promising. We expect the 
shape of its power spectrum to be independent of Or in 
an integrable system, but not in a chaotic system. 

We thank M. Brack, A. V. Khaetskii, and M. Plc- 
tyukhov for stimulating discussions and P. Brouwer for 
a helpful clarification. The work has been supported by 
the Deutsche Forschungsgemeinschaft (OZ and KR) and 
the EU Spintronics Research Training Network (DF). 



Our numerical simulations show that in integrable sys- 
tems both the spin direction and the phase oscillate al- 
most periodically during the orbital motion. Therefore, 
after a short transient period, M.(L) usually saturates. 
One exception we found is the circular billiard. Here, ow- 
ing to angular momentum conservation, all trajectories 
efficiently accumulate area, and M.(L) ~ 2sm(x)/x for 
Or <C 1, where x = 8^Lr/2L^ and r is the radius [f3 |. 

Figure [3 shows the relative quantum correction to re- 
flection, sn/snw, as a function of Or for chaotic and 
integrable geometries (at B = 0). Positive (negative) val- 
ues of SIZ/SIZ^ ) indicate WL (AL). For chaotic systems 

SR/SRM = (5K diag + SK loop )/(5K2 g +Sn[ J op ) = (M~ v )l 
is given by Eq. ©■ For the numerical calculation of 
A4(L) only (backscattered) orbits starting and ending at 
one lead are considered (since they are closed, the ini- 
tial spin relaxation is reduced compared to non-closed 
paths). The chaotic DS (solid curve) and DD (dashed 
curve) billiards show a very similar WL-AL transition 
with increasing Rashba strength. The escape length in 
units of Lb is P c / (w+w'), where P c is the perimeter of the 
cavity. Hence, given Lr, one can also conclude that AL 
is absent in smaller quantum dots [for fixed P c /(w + w') 
or w + it)'], as supported by experiment 0. 

The results for the integrable QC billiard (dashed curve 
with circles) are based on a numerically obtained length 
distribution, which is no longer exponential. The tran- 
sition to AL in the integrable billiard is much less pro- 
nounced and occurs at clearly higher Or, compared to its 
chaotic counterparts owing to the slower spin relaxation. 
Hence there exists an extended regime of SO strengths, 
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